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We consider Yukawa theory in which the fermion mass is induced by a Higgs like scalar. In 
• our model the fermion mass exhibits a temporal dependence, which naturally occurs in the early 

I Universe setting. Assuming that the complex fermion mass changes as a tanh-kink, we construct an 

exact, helicity conserving, CP-violating solution for the positive and negative frequency fermionic 
mode functions, which is valid both in the case of weak and strong CP violation. Using this solution 
we then study the fermionic currents both in the initial vacuum and finite density /temperature 
setting. Our result shows that, due to a potentially large state squeezing, fermionic currents can 
exhibit a large oscillatory magnification. Having in mind applications to electroweak baryogenesis, 
we then compare our exact results with those obtained in a gradient approximation. Even though 
the gradient approximation does not capture the oscillatory effects of squeezing, it describes quite 
well the averaged current, obtained by performing a mode sum. Our main conclusion is: while 
the agreement with the semiclassical force is quite good in the thick wall regime, the difference is 
sufficiently significant to motivate a more detailed quantitative study of baryogenesis sources in the 
^ . thin wall regime in more realistic settings. 

PACS numbers: 98.80.-k, 04.62.- 
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I. INTRODUCTION 



Electroweak baryogenesis [l| is a very appealing idea, and yet the mechanism for dynamical baryon creation at the 
electroweak scale has suffered some serious blows. Firstly, in the mid 90s it was found that the electroweak phase 
. transition in the standard model is a crossover While at first supersymmetric extensions looked promising, the 

most popular supersymmetric model - the MSSM - is almost ruled out on two grounds (a) it cannot give a strong 
enough phase transition for the observed Higgs mass ^ and (b) it cannot produce enough baryons consistent with 
electric dipole moment bounds fjl-flpl ( albeit in some models resonance between fermionic flavors can be helpful 
I to increase baryon production Jw, [l3l - [l5| ). The models that are still viable are the supers ym metric models with 
5—1 • additional Higgs singlet (s) [l6l.ll7l| both because they allow for a stronger phase transition [l8l420l | and generate more 
baryons j2l| - |25| . In addition, general two Higgs doublet models [H, [13 and composite Higgs models j28l - l30| are still 
viable. Works on cold electroweak baryogenesis [3ll - l33| are also worth mentioning. In summary, while electroweak 
baryogenesis has been a very attractive proposal, precisely because it is testable by contemporary accelerators, recent 
experiments have cornered it to models where most researchers have not focused their attention during the pre-LHC 
era. Hence, at this stage theoretical work, that will refine our ability to make a quantitative assessment of electroweak 
baryogenesis in different models, is still a worthy pursuit. 

One of the most important unsolved problems in dynamical modeling of electroweak baryogenesis is a reliable calcu- 
lation of the CP- violating sources that bias sphaleron transitions [13, , which at high temperatures violate baryon 
number. In the fermionic sector the most prominent CP-violating source is the fermionic axial vector current [36l.l37|. 
since that current directly couples to sphalerons, and can thus bias baryon production. There are essentially two 
approximations used in literature to calculate axial vector currents: 



(a) the quantum-mechanical refiection [38|-|40| used in the thin wall case, and 
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(b) the scmiclassical force [O, [H, HI EHill used in the thick wall case. 

In general thin wall baryogenesis is more efficient in producing baryons. Its main drawback is that the calculational 
methods used are unreliable: one calculates the CP-violating reflected current ignoring the plasma, and then inserts it 
into a transport equation in an intuitive (but otherwise rather arbitrary) manner [40| . How bad the situation can get 
is witnessed by the controversy that developed around the work of Farrar and Shaposhikov [4^ (who used a quantum 
mechanical reflection to calculate the source). The subsequent works (45l - [47| came up with an orders of magnitude 
smaller answer for baryon production. And yet these latter works used unreliable methods that e.g. violate unitarity, 
such that the issue remained unsettled ^. So, the problem of the source calculation in the thin wall regime remains 
still to a large extent open. 

In the thick wall case the situation became much more satisfactory after the works of Joyce, Kainulainen, Prokopec, 
Schmidt and Weinstock [i^ IH, . It was shown that one can calculate the scmiclassical force (which rather straight- 
forwardly sources the axial vector current) from first principles and in a controlled approximation from the Kadanoff- 
Baym (KB) equations for Wightman functions. These KB equations are the quantum field theoretic generalization 
of kinetic equations. The positive and negative frequency Wightman functions represent the quantum field theoretic 
generalization of the Boltzmann distribution function, that provide statistical information on both on-shell and off- 
shell phase space flow. In a certain limit, when integrated over energies, the Wightman functions yield Boltzmann's 
distribution function. When written in a gradient approximation, the KB equations can be split into the constraint 
equations (CE) and the kinetic equations (KE). The authors of Refs. [i^, have rigorously shown that, in the 
presence of a moving planar interface, in which fcrmions acquire a mass that depends on one spatial coordinate, single 
fermions live on a shifted energy shell, which to first order in gradients (linear in h) and in the wall frame equals 



where m(z) = mji{z) -\- imj (z) = |m|(z)e'^(^) is the fer mion mass, which varies in the z-direction in which the wall 
moves, k is particle's momentum, k± is the momentum orthogonal to the wall, and s = ±1 is the corresponding spin. 
This energy shift acts as a pseudo-gauge field (also known from condensed matter studies), which lowers or increases 
particle's energy. Relation ([Ij clearly shows that particles with a positive spin orthogonal to the wall and a positive 
frequency (as well as particles with a negative spin and a negative frequency) will feel a scmiclassical force that is 
proportional to the gradient of 6* = arg[m]. Particles with a negative spin and a positive frequency will feel an opposite 
force. This force appears in the kinetic equation for the Boltzmann-like distribution functions f±s, and reads 

F,. = -M!i£±..^iMML. (2) 

It was also shown that this force then sources an axial vector current, which in turn can bias sphalerons. 

The work of [l^ l47l . Isol - fs^ has shown that, in the case that fcrmions mix through a mass matrix, there is an 
additional CP-violating source resulting from flavor mixing. This was put on a more formal ground by [l^, where 
a flavor independent formalism was developed, and where it was shown that flavor non-diagonal source is subject to 
flavor oscillations induced by a commutator term of the form i[M, /], not unlike the famous flavor (vacuum) oscillations 
of neutrinos. This idea was further developed by [s^. Since we do not deal with flavor mixing in this work, we shall 
not further dwell on this mechanism, which should not diminish its importance. In passing we just mention that in 
most of the relevant parameter space of e.g. the chargino mediated baryogenesis in the MSSM, the scmiclassical force 
induces the dominant CP-violating source current 

We shall now present a qualitative argument which suggests that in many situations thin wall sources can dominate 
over the thick wall sources (calculated in a gradient approximation). If true this means that any serious attempt to 
make a quantitative assessment of baryon production cannot neglect the thin wall contribution. To see why this is so, 
recall that a gradient approximation applies for those plasma excitations whose orthogonal momentum, k± = 27r/A_L, 
satisfies: 

27r 

> — (THICK WALL) , (3) 
L 



Research on the topic subsided not because the problem was resolved, but because standard model baryogenesis was ruled out based on 
equilibrium considerations alone Q|. 
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where L is the typical thickness of the bubble wall. On the other hand, the thin wall approximation belongs to the 
realm of momenta which satisfy 

277 

A:^ < — (THIN WALL) . (4) 

ij 

Typical momenta of particles in a plasma (per direction) is k±_ ^ T. Now, unless LT ^ 2tt, we have a larger 
or comparable number of particles in the thin and thick wall regimes! But, since the thin wall source is typically 
stronger, unless thermal scattering significantly suppress the thin wall source, it will dominate over the thick wall 
source. It is often incorrectly stated in literature that the number of particles to which thin wall calculation applies 
is largely phase space suppressed, i.e. that their number is small when compared to the number of particles to which 
the semiclassical treatment applies. So, to conclude, it is of essential importance to get the thin wall source right if 
we are to claim that we can reliably calculate baryon production at the electroweak transition in a model. 

We believe that this represents a good motivation for what follows: a complete analytic treatment of fermion tree- 
level dynamics for a time dependent mass. The time dependence has been chosen such to correspond to a tanh-kink 
wall, because it is known that this represents a good approximation to a realistic bubble wall |M . [ssj , and equally 
importantly, in this case one can construct exact solutions for mode functions. Before we begin our quantitative 
analysis, we recall that a related study for the CP even case and planar wall has been conducted by Ayala, Jalilian- 
Marian, McLerran and Vischer [56j . while a semianalytic, perturbative treatment of the CP- violating case has been 
conducted in Ref . [s^l • The main advantage of the latter study is that it allows for a general profile of the CP-violating 
mass parameter, the drawbacks are that the method is semianalytic (the final expression for the source is in terms of 
an integral), and furthermore it is perturbative, such that it can be applied to small CP violation only. To conclude, an 
exact treatment of fermion dynamics in the presence of a strong CP violation is highly desirable, and this is precisely 
what we do in this paper. 



II. THE MODEL 



Here we consider the free fermionic lagrangian of the form, 

Co = ■ipi'y^d^ip - m*{ijRi}L - rmpLipR , (5) 

where ipL = Pl'^P and i/jr = PriP are the left and right handed single fermionic fields, Pl = (1 — 7^)/2 and Pr = 
(1 + 7^)/2 are the left and right handed projectors, and 7^ and 7^ are the Dirac gamma matrices. We shall assume 
that the fermion mass m is complex and space-time dependent. This can be generated e.g. when a Yukawa interaction 
term, Cy = —y4>ipLipR + h.c, is approximated by —y{(l))ipLipR + ^-c, where (</>) stands for a Higgs-like scalar field 
condensate which can generate a space-time dependent fermion mass, 

■m{x) y{(f){x)) , (6) 

where y is a (complex) Yukawa coupling. The Dirac equation implied by ([5]) is 

ij^dfj_ip — 'm*'ipL — niipR^ = . (7) 

In this paper we consider the simplest case: a single fermion in a time dependent, but spatially homogeneous, 
background. Such situations can occur, for example in expanding cosmological backgrounds [58[, or durin g se cond 
order phase transitions and crossover transitions in the early Universe. In this case helicity is conserved [59l - l6l| . We 
shall perform the usual canonical quantization procedure, according to which the spinor operator ip{x) satisfies the 
following anti-commutator {h = 1), 

{lP^{x,t),^pl{x',t)} = S^pS^x-x') . (8) 

In the free case under consideration, the Dirac equation (jT)) is linear, and consequently ipi^) can be expanded in terms 
of the creation and annihilation operators, which in the helicity basis reads, 

^(f , t)= [ y'''HhXh{l t) + e-^^^-%\v^{k, t)] , (9) 

•' ^'^^> h=± 

where Xhik,t) and Vh{k,t) are particle and antiparticle four-spinors. a^^ and bj:^^ are the annihilation operators that 
destroy the fermionic vacuum state |f2), a^^jri) = = fo^r^jri), while at^ and bt^ are the creation operators that create 
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a particle and an antiparticle with momentum k and helicity h. These operators obey the following anticommutator 
algebra, 



where all mixed anticommutators are zero. The momentum space quantization conditions (jlOp and the position space 
quantization rule ([8]) have to be mutually consistent. This imposes the following consistency condition on the positive 
and negative frequency spinors, 

^[Xha{k,t)xlp(k,t) + Vha{-k,t)vlp{-k,t)] = 5ap ■ (11) 

h=± 

This is usually supplied by the mode orthogonality conditions, 

Xh(k, t) ■ Vh{k, t) = Q = Vh{k, t) ■ xh{k, t) . (12) 
and by the mode normalization conditions, 

xl{k,t) ■ Xh{k.t) = 1 = ul{k,t)-Vh{k,t), (13) 

which - as we will see below - are chosen to be consistent with the more general requirement (|lip . Although the 
orthogonality condition ([T2l) is usually met, it is however not a necessity. Important is that the mode functions span all 
of the Hilbert space, which is true in this case. Because we consider a system which is time-translationally invariant, 
helicity is conserved, and it is thus convenient to work with helicity conserving spinors 

^ ( *:t'i! ) ^ ( fiJiS! ) <") 

where £,h{k) is the helicity two eigen-spinor, satisfying h^h = h£^h, where h = k ■ a is the helicity operator and /i = ±1 
are its eigenvalues. 

We shall work here with the Dirac matrices in the chiral representation, in which 



where the last equalities follow from the usual direct product (Bloch) representation of the Dirac matrices. Here 
and cr* are the Pauli matrices obeying, p^p^ = (5*^ + ie'-^'p' and crV-' = (5*^ + ze*^'cr'. The left and right projectors are 
then, 



1-7 f I 0\ 1+P r n 1+7' fO 0\ 1-P . 



which can be used to write, ipL = PliP and ipfj, = Pb.iP as it is done in ([5]-[7]). Now, making use of Eqs. (pi-[T5)) in the 
Dirac equation ([7]) one gets the following four equations for the component functions 

iLh + hkLh = mRfi 

iRh — hkRh ~ m* Lh (17) 

and 

iLfi — hkLh ~ "niRh 

iRh + hkRh = Tn*Lh , (18) 
where the mass can be complex and time dependent, m ~ Tn{t), and the modes are normalized to unity, 

\Lh\^ + \Rh\' = l = \Lh\' + \Rh\^ (19) 
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The equations of motion for and Rh can be decoupled, resulting in the second order equations, 

Lh + uj'^Lh - —{Lfi - ihkLh) = 



771 



TO 



Rh+u'Rh -{Rh + ihkRh) = 0, (20) 

TO* 

where = fc^ + |TO(i)p. For the case at hand a better way of proceeding is to go to the positive and negative 
frequency basis, defined by: 

u±h = ^{Lh±Rh). v±h = ^{Lh±Rh), (21) 

since then the equation of motion can be reduced to the Gauss' hypergeometric equation. Indeed, from ([T7l - [T8|) 
and ([211) it follows, 

iu±h T mjf{t)u±h = -{hk±tmj)u^h 

iv±hTmR{t)v±h = (hk^imi)v^h, (22) 
which, when decoupled, yields a second order equation, 

•■ ^ ■ , /^;2 , I |2 , • , ruR^i \ n /oq^ 

u±hT^-ri~> u±h+\k + |?7i| ±imR + -—— \u±h^'J- (23) 

hk ± mil \ nk ± imj I 

So far our analysis has been general, in the sense that we have assumed no special time dependence in 'm[t). In order 
to make progress however, we have to make a special choice for m(t), which is what we do next. 



III. MODE FUNCTIONS FOR THE KINK PROFILE 

In Ref. [56j an exact solution of the Dirac equation was found for a wall of arbitrary thickness with a kink wall 
profile cx tanh(— z/L), where L = l/\ characterizes the wall thickness. Here we generalize this solution to include 
CP violation. While in this paper we consider only a time dependent mass profile, the generalization to the planar 
(^-dependent) case is straightforward, and will be considered separately. Constructing an exact solution is important 
for baryogenesis since one can then consider in detail how the CP-odd quantities that source baryogenesis (directly 
or indirectly) depend on the mass profile, and in particular investigate what is the optimal profile and its duration. 
Unfortunately, analytic solutions cannot include plasma scattering and width effects, whose treatment will be therefore 
typically left to numerical simulations. 

Here we assume the following 'wall' profile 

m{t) = nil + fTi'2 tanh ^ ^ , (24) 

where r = I/7 represents the time scale over which the wall varies (for convenience we shall use the terms 'wall' and 
'profile' interchangeably). Both toi and TO2 are complex mass parameters. In the case when a single Higgs field is 
responsible for the phase transition, one expects that both real and imaginary part of m(t) exhibit a similar behaviour, 
which is reflected in the Ansatz ([24]) . Moreover, we do not know how to construct exact solutions when different time 
scales govern the rate of change of the real and imaginary masses. Nevertheless, we believe that the Ansatz ([24]) 
represents quite well realistic walls for a wide variety of single stage phase transitions, c/. Refs. [s^. [55|. 

Note that the thin wall limit is r — > (7 — > 00). In that limit the mass function becomes the step function 
Ansatz (jBip . whereby m± = toi =p 1712- In appendix El we construct the normalized fundamental solutions of Eqs. 
P7p for a constant mass. The thin wall case is treated explicitly in appendix |B1 The thin wall results serve as a check 
for the kink wall case in the appropriate limits. Moreover it allows for a quantitative comparison of the thick wall to 
the thin wall results. 

Since the ratio of the real and imaginary parts of the mass mi{t)/mR{t) is time dependent, the Ansatz (|24p contains 
CP violation (which can be either small or large, depending on how much the ratio mi{t) /muit) changes. Since the 
physical CP-violating phase is in the relative phase between mi and TO2, one can perform a global rotation of the left- 
and right-handed spinors that does not affect CP violation. It turns out that the equations of motion simplify if one 
performs a global rotation that removes the imaginary part of 7712. The constant rotation that docs that is 

m{t) m{t)e''^ , x = arctan ( -IUlL | . (25) 

V m2Bj 
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In that case 

mi = mm + mu , 1122 = m2R . 
This rotation is important, because the mode equations (|23p significantly simphfy to become 

u±h + (t^^(i) ± imR)u±h , 



(26) 



where uj'^{t) = fc^ + |m(t)p. Furthermore, from p2|) one can infer that v±h obey the same equations as u±^h- In 
what follows, we show that these equations can be reduced to the Gauss' hypergeometric equation. 



To show this, it is instructive to introduce a new variable. 



1 1 



- — - tanh ( 



(27) 



in terms of which 

m{t) ^ mi + m2(l — 2z) , 

with 7 = 1/r. Eq. ((26|) becomes, 

47^[z(l-z)]2^+472(l-2z)z(l-z)^ 



mR{t) = -2m2RZ = - 



7TO2_R 



cosh^(-t/T) 



'A'jm2Rz{l - z) 



P + m] + {miR+m2Rf - AzmiRm2R - 4z(l -z)m2fl(m2i? ± 27) 



u±h = . 



Now, performing a rescaling. 



and choosing 



where 



u±h = z'^{l^zfx±h{z) 



'2T' 



'2~ 



Wzp ^ tj(t =Foo) = y fc2 + + (mi/j ± m2i?)2 , 
yields the following Gauss' hypergeometric equation for x±hi 



z{l - ^)-^ + [c - (a± + b± + 1)^]^ - a±^± 



where 



a±=a + /3 + l=F« 



m2_R 



&± = a + /3 ± I 



m2fi 



X±h{z) = 0, 



c = 2a + 1 , 



(28) 



(29) 



(30) 



(31) 



(32) 



(33) 



Note that the rescaling (|29p was chosen such to remove the terms oc l/zandoc 1/(1 — z) from Eq. ([5^ . Since a±,6±,c 
are non-integer, the two independent solutions for x±h are the usual ones. A detailed normalization procedure is 
provided in Appendix [D] and the result arc the following normalized early time mode functions 



U+h 



(1) oj- + {miR + m2R) an \0 IP r 1. ^ 
= u\l^ = J ^— X z"(l - zf X 2Fi{a+,b+;c;z) 



(1) 

u-h = u'_j^ = . X 



hk — imi 



- {miR + m2R) 
2a;_ 



z"(l-z)/5 X2Fi(0 



(34) 
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These functions are valid of course for all times. They are called early time mode functions because at early times 
{t — >■ — oo) they reduce to the positive frequency mode functions (|D2p . and they are normalized as, + = 1, 

which follows from Eqs. dH]) and ([HI), see also Eq. (jPlSp . 

For completeness, we also quote the second pair (jDl[) of early time solutions, 

p+c-a+-b+ ^ ^p^^^ - a+, 1 - b+- 2-c-z) 



(2) ^ hk-imi ^ /wT+T^J^Tn^ i_ 



X2Fi{l-a-,l-b-;2-c;z). (35) 



Just as before, at early times {t —oo) these solutions reduce to the negative frequency mode functions (jD2l) . and 
they are also normalized as, jit^^p + [m'^^P = 1. 

An analogous procedure as above yields the following normalized fundamental solutions suitable for late times, 



y 2w+ 



u^ll = -^^?=^ X >+ - - X ."+^-(1 - .)^+-— X ,F,(1 - a_, 1 - 2 - 5; 1 - z) (36) 



and 



- ■' X z°'{l - z)P X 2Fi{a+,b+;c;l - z) 



u 



-h 



2uj+ 

.(2) hk-imi UJ+ + {mm - m2R) , , ~ -, n /o-tn 

u'_i = — === X J X z°'{l- zf X2Fi{a-,b_;c;l- z), (37) 



while the late time solutions ([57]) reduce at asymptotically late times to positive and negative frequency solutions 
cx eT*"+*, respectively, see Eq. (jD7|) . 

Now, a general early time solution can be written as a linear combination of the fundamental solutions p4H35p : 
for simplicity we shall take here p4p for the early time solutions. Similarly, general late time solutions are a linear 
combination of the fundamental late time solutions 



u±h = a±/iW±^ + P±hU'±h ' (38) 

where a±h and P±h arc complex functions of k (for spatially homogeneous systems they arc functions of the magnitude 
only) that satisfy the standard normalization condition, 

\a±h? + \P±h? = 1- (39) 

Now, upon choosing p4p as the early time solutions and making use of the matching between the general early and 
late time solutions 

u±h{k,t) =u±h{k,t) (40) 
and of the relation for the Gauss' hypergeometric functions (jPSp one gets, 

a±h = 



/ UJ+ [uj^ ± (mi + m2fl)] r(c)r(a± + b± - c) 
uj- ± (mm - m2fl)] r(a±)r(6±) 



^ _^ ± {'miR + m2R)] r(c)r(c - a± - b±) ^^^^ 

y LO-[uj+ =F {miR - m2R)] T{c - a±)r(c - b±) 



It can be shown that 



a+h = a-h 



(42) 
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Useful identities here are 



uj- T {miR + m2R) = 



±ijj\ T T 2m2ii)^ 
=Fw^ ± (w+ ± 2m2fi)^ 



4m 



2i? 



(43) 



Because a±;i and /3±;i are functions of a±, 6± and c, (just as in the thin waU case (jB6IIB7|) ) tliere are no CP 
odd contributions in the mode mixing (Bogohubov) coefficients (|4ip . a±h and I3±h are indeed the usual Bogoliubov 
coefficients that transform an asymptotically early time vacuum state to a late time vacuum state. Hence n±h = |/3±h P 
is the particle number observed by a late time observer, in the late time state that evolves from the early time positive 
frequency vacuum state. 

To make contact with the thin wall case (|B6p . we take the limit 7 oo in (|4T|) to get, 



I [uj^ ± [miR + m2R)] 



T m2R 



(44) 



It can be checked that Eq. (|44p satisfies = Moreover, since a;„ — aj+ — 2m2R < 0, the and (3-fi are 

always positive. One can show that a±h and (3±fi given in (|4ip obey |a±/jp + |/5±/tP = 1, as they should. This equality 
follows from. 



\a±h\ 



n±h = \l3±h\ 



sinh sinh 
sinh ( £[i£^^f£±±H!22Hl'| ginh ( "'"+-"-+''"^ 



27 



27 



sinh sinh 



(45) 



from which it also follows that |Q:±?ip — 1 ~ |/3±/tP- Now, taking a thin wall limit 7 — >■ 00 in (P5|) yields 



m_|_p — (w- — w+)^ 



(46) 



where we made use of Given that ^nri^j^ = \m- — m+p. This expression agrees with the thin wall particle number 
Eq. (jB7|) derived in appendix [Bl 

It is interesting to note that, although particle number agrees, the Bogoliubov coefficient (3±h in the thin wall limit 
(HH) appears very different from the one derived explicitly for the thin wall (|B6p . For instance, the coefficients in (|B6p 
are complex and depend explicitly on hclicity, whereas the limiting coefficient (j44p is real and helicity independent. 
A similar situation occurs for a±;t, sec (|E8p . The apparent discrepancy is caused by an overall phase factor by which 
the coefficients in the thin wall limit differ from those directly computed for the thin wall. This phase factor does not 
affect particle number and can be removed by a global rotation of the (anti)particle spinors. In appendix lEl we show 
explicitly how the kink wall case and thin wall case are connected. 

The particle production can also be analyzed in the opposite limit, 7 — > 0. In this thick wall regime particle 
production is exponentially suppressed as. 



7-!-0 

n±h — > exp 



2m2fl) 



7 



(47) 



which is also what one expects. However, note that when 7r(a;+ + w_ — 2m2fl) < 7, the suppression is not large. This 
is demonstrated in figures [T] and where the particle number is shown as a function of k for several different wall 
thicknesses. In figure[T]the mass parameters are toik = m2R and to/ <C toi/j, m2R. In this case CP violation is weak. 
For these mass parameters the thin wall particle number (|46|) . represented by the dashed line, reaches the maximal 



particle number n±h = 2 as fc — )■ 0. For thicker walls (decreasing 7) the particle number is exponentially suppressed 
with respect to the thin wall. For very small k the suppression is much smaller, since t^{ijJ+ + — 2TO2fl)/7 ^ 

In figure [2] the mass parameters are chosen such that mj,miii <^ m2R. In this case CP-violation is maximal for 
the thin wall in the limit fc — >■ 0, see also (|B9p . The maximal particle number in this limit is 1, which indicates an 
inverse population. This inverse population, induced by large CP violation, is a novel result and, as far as we know. 
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0.5 r 




k 



2 4 6 8 10 



Figure 1: Plot of n±h ~ \l3±hf function of k in units 
of m2R- Parameters are tuir — m2R, mi = Q.lm2R. 
The dashed line is the thin wall solution of n±h . The 
other (solid) lines show - from top to bottom - the particle 
numbers for 7 — 10m2R (blue, dark), 7 = m2R (red) and 
7 = 0.1m2R (orange, light). In general particle number 
is suppressed for decreasing y (a thicker wall), but still a 
large particle number is reached when k, mi <C m2R and 
miR — m2R ~ 0. 



n±i, 




2 4 6 8 10 



Figure 2: Plot of n±h ~ |/3±hP as function of k in units 
of m2R. Parameters are mi = 0.1m2R, miR — 0.1m2R. 
The dashed line is the thin wall solution of n±h.. The 
other (solid) lines show - from top to bottom - the particle 
numbers for 7 — 10m2R (blue), 7 = m2R (red) and 7 — 
0.1m2R (orange). Because miR — m2R < 0, an inverse 
population is reached. The CP-violating phase is maximal 
because m-^- ~ — m_ . 



not noticed in literature before. For thicker walls the particle number is still suppressed, but much less than for the 
mass parameters in figure [1] In fact, for mm = mj = the particle number is unsuppressed in the limit fc ^> 0. 

A large late time Bogoliubov particle number for a free fermionic system indicates large squeezing. It is interesting 
to see what effect such a large squeezing may have on the fermionic currents. In particular, we are interested in the 
CP-odd fermionic axial vector current that couples to sphalerons. The next section is devoted to computing these 
currents in the setting of a tanh-kink wall. 



IV. THE CURRENTS AND CP VIOLATION 

In this section we consider the evolution of the two point Wightman functions, defined as the expectation values [H, 

Hi 

^-^Q/s" = l-^afsi'^^'") = -('/'/3(w)'0q(w)) ; lS~^{u,v) = lS>p{u,v) = {lpa{v)i) fj{v)) , (48) 

and which satisfy the homogeneous Dirac equations ([7]) 

(i7^a^ -mn- imn^)iS^J (u, v) = G . (49) 
For the problem at hand, when written in a Wigner mixed representation 

v)= j ^o^'-^^-^hS^^ik; x), (x = {u + v)/2) , (50) 
the fermionic Wightman function can be written in a helicity block-diagonal form 

iS+-{x-k) = iS< ^ ^S<, ~iy'>S+- ^{p''gah)<^\{l + hk-a), (51) 

where a"^ , p'^ [a ~ 0,1,2,3) are the Pauli matrices and gah are the (off-shell) distribution functions measuring the 
vector, scalar, pseudo-scalar and pseudo- vector phase space densities of fermions, respectively. Their on-shell version 

/J 7 
-^gah-, (a = 0,1, 2, 3) (52) 
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satisfy the following equations of motion [36|, |59| , 



foh = 

fih + 2hkf2h - ^mif^h = 

/2/1 - 2hkfih + ImfiJvi = 

/3/1 + 2mifih - 2miif2h = , 



(53) 



where here k = \\k\\. To make the connection with section Hill and Appendix [Bl we note that one can express fah in 
terms of u±ii or Lh and Rh as follows ^: 

foh = \u+h\'' + \u^h\^ = \Rh\^ + \Lh\^ ; hh = 2^[u+hu\] = \Lh\^ - \Rh\'' 

fih = \u^h\^-\u+h\^ = -2^[LhRl]; f2h = 2^[u+hu*_^] = -2'^[LhRl] . (54) 

such that fih + 1/2/1 = —2LhR\- From Eq. ()A9P and ([54| we immediately obtain that for t —00 {z 0), 



/o/i ~ -'^ '■ 



_ _hk 
hh - ^_ 



where we took account of u+hU*_i^ = —{kh + iTO/)/(2w_). z = exp(2i/r)/[l + exp(2t/r)] exp(2t/T) (as t 
and of 2^1 (a, b; c; 0) = 1. Inserting Eqs. (|55p into the particle number definition [59|, 



nh{k,t) 



mnfih + m//2/i + hkf^h 1 
2w 2 



(55) 
—00) 

(56) 



yields that nh{k,t) = for t — >■ —00, as it should be since we have prepared the initial state to be in the pure free 
vacuum. 



One can also consider the statistical particle number [62[, 

nh± = ^foh ± \ \l flh + flh + fi 



ih ■ 



(57) 



A statistical particle number is defined as the particle number associated with the basis in which the density operator 
is diagonal j62| . Statistical particle numbers can be used as a quantitative measure of state impurity, i.e. of how 
much a state deviates from a pure state. From previous work we have learned that the statistical particle number is 
constant in the absence of interactions. This can also be seen from the kinetic equations (|53p . which give 



j^{fih + .fih + fk)=o- 



(58) 

Of course, when interactions are included, the righthand side of Eqs. (|53p is in general nonzero. Here we consider 
a free Dirac equation (|49p . and therefore the statistical particle number should remain constant. Indeed, Eqs. ((55)) 
imply 

.f?h + fk + fiH=[\u+H\' + \U-h\']' = l, 

such that the statistical particle numbers of a pure state are trivial. 



foh ± \[?L + llh^~7ih 

Thus the statistical particle number is either or 1, the latter corresponding to a fully occupied Dirac sea. 



Note that, duo to diflfcroncc in conventions, there are sign differences when compared with Ref. |59|| . 
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The exact solutions for the phase space densities fah as in Eq. ((54)) are comphcated functions containing products 
of hypergeometric functions, which can be analyzed numerically. Analytically, we can study the behavior of fahS in 
certain asymptotic limits. By making use of Eqs. (p6l - [4T|) . we get at asymptotically late times. 



~ >-oo 

U+h — > Ct+h 



I u}+ + {miji - m2R) 



2w4 



I UJ+ - {mm - m2fl) 



u-h 



-a-h 



hk - mil / w+ - [mm - m2R) _ 



VP 



2w_, 



hk — tmj / cj-i- + {miu — m2R) 



(59) 



From these and Eqs. (|39|) . p2|) and ([54)) we easily obtain, 
1 



/o/i 
f+ 

Jlh 
fth 



rriR 
mj 

2hk sin(2w+t) 



(i-2i/3±n 



2^yk'^ + mj 
2 cos(2aj+t) 



l3h 



y^k"^ + mj 
— (l-2|/3±n 

UJ+ 



W-l_-\/fc^ + mj 
2 cos(2cj+i) 



[5R[«±,./3JJ cos(2t^+t) + ^[a±hP±h] sin(2L^+t)] 
[^[a±hf3'^h\'^i^[m+] + ^[a±hf3^,^]hkuj+] 



[^[a±hf3±f,]^ik^[m+] - '^[a±hl3^,^]uj+mi] 



Ld+y/k'' + mj 



(60) 



Since we are primarily interested in the CP-violating axial vector currents, which can bias sphalerons, here we shall 
focus our attention primarily to fsh ■ Wc can compute the CP-odd and CP-even axial vector phase space densities as 
follows. 



h=± 



4|a±^||/3;jm/ 

, sm(2w+< - Aip) 

yjk^ + mj 



-(1 - 2|/3±| ) + — cos(2w+t - A(^) 



where 



A(/3 = - (p;3 



(61) 



(62) 



(^Q, and ipp, are the phases for a±h and /3±;i in Eq. (jH)). In the thin wall limit 7 cx) the phases of oi±h and are 
zero (see Eq. ()E8P ) and the CP-odd and -even phase space densities coincide with those in Eqs. ([Clip (in the free 
vacuum). In the opposite 'thick wall' limit, 7^-0, 



with 



VP 



>o 1 



7 

w_ + UJ. 



a±h — > e 
/3±/i — exp 



UJ- log(w_) — LOj^ log(cL'_|_) 



7r(aj-|_ + w_ — 2m2R) 
2^ 



(63) 



>o 1 



7 

w_ — UJ. 



m2fi j log 
w_ log(w_) + w+ log(w+) 

"^2i? log 



"12_R log 



+ m-2_R lo. 



m2fl 



(64) 
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Figure 3: Odd (solid blue, darker) and even (solid 
red, lighter) part of the exact solution for f^h for 7 = 
lGm2R. Parameters: k = 0.1m2R, 7 = 10m2R, mj = 
0.1m2jj, miR = 0.17712/?. The thin wall solutions (dashed) 
are constant for t < 0, and for t > they oscillate with a 
frequency 2oj+. The thin wall and exact result are nearly 
identical. 



Figure 4- Odd (solid blue, dark) and even (solid red, 
lighter) part of the exact solution for f^h for 7 = 
m2R. Parameters: k = 0.lm2R, 7 ~ m2R, mi = 
Q.lm2R, miR — 0.1m2R. For this wall of thickness 
7 = rn2R, the amplitude of the even and odd part of fsh is 
slightly larger than the thin wall case (shown as dashed). 
Also, there is a moderate phase shift compared to the thin 
wall. 



Thus, in this hmit the phases ipa and —ip/j grow hncarly with I/7. In general, the oscihating CP-odd and CP-even 
phase space densities in the thick wall case therefore experience a phase shift compared to the thin wall, which can be 
large for small 7. Moreover, their amplitudes, which are proportional to |/3±hP (sec Eq. (|61l) . are generally suppressed 
compared to the thin wall limit. However, in the previous section we have demonstrated that the amplitude |/3±/i|, or 
particle number nh± , can become much less suppressed for a certain choice of parameters that leads to large squeezing, 
see figures [THU Due to this state large squeezing, the oscillations of the phase space densities are amplified. 

Examples of this enhanced oscillatory behavior are depicted in figures [3HS1 Here we show the exact solution for 
the odd and even part of fsh using Eqs. ([54]) with solutions (p4)) . and compare to the thin wall solutions (jClip . The 
parameters are chosen such to generate large squeezing, thus k ^ m2R and the mass parameters are the same as 
those in figure O Close to the thin wall limit (7 = 10m2R, figure [3]) the exact solution for f^h almost coincides with 
the thin wall result. For a thicker wall (7 = m2R, figure U]) the amplitude of oscillations remains large and there is a 
modest phase shift compared to the thin wall. Finally, for a thick wall (7 = .lm2R, figure[5]) there is a large oscillatory 
enhancement and a large phase shift. 

For thick walls a large squeezing can thus give a large enhancement of oscillations of the phase space density /a^. 
However, due to the fact that phases of different modes k differ, the oscillatory behavior may (partially) disappear 
when the corresponding currents are computed. These currents are related to the phase space densities fah as 

/(Pk 
j^^fah{k;x) , (65) 

where jo/i and j^h denote the 0th components of the vector and axial vector current, respectively, and jih and j2h 
arc the scalar and pseudoscalar densities, respectively. In the following section we compare these integrated currents 
to those computed in a gradient approximation. First however, we shall show how to compute the currents for more 
general initial states. 



A. Generalized initial state 



So far the initial state has been taken in the free vacuum, such that nh{k, t) = Q for t — > —00. We can also consider 
the initial state to be thermally occupied, such that 

nh{k,t) *^~°°> nth = fl^'^ , 1 • (66) 
eP"- + 1 

The initial phase space densities that give this initial thermal state via Eq. ([55)1 are now 

/o. = l; /r. = -^^(l-2nth); /-,--^i^(l-2nth); /3,^- — (l-2nth), (67) 

UJ— UJ— UJ— 
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Figure 5: Odd (solid blue, dark) and even (solid red, 
lighter) part of the exact solution for f^h for 7 = 
0-lm2R (solid line). Parameters: k — 0.1m2R, 7 = 
0.1m2R, mi = 0.1m2ij, mm — Q.lm2R. When CP- 
violation is maximal and squeezing is large (see figure 
there is a large oscillatory enhancement for the odd 
and even phase space densities m the thick wall regime 
(7 = G.lm2R). Moreover, there is a large phase shift 
compared to the thin wall result (dashed). 



Moreover, the statistical particle number ([57]) for t —00 is 

nh+ = 1 - "th 
fih- nth ■ 

Using these relations the currents can also be written in terms of statistical particle number 



(68) 
(69) 



flh = {^hA 



nil-); f: 



2h 



-{nh+~nh-); f. 



3h 



hk 
[rih- 



rih 



(70) 

Here, one can also consider general initial densities. In that case the statistical particle numbers nh+ and rih- are 
both in the range [0, 1] and conserved, but they are not necessarily related as in an initial state that is thermally 
occupied. The scalar density foh is in the range [0, 2] and is still conserved. The initial state may be called an initial 
dense state, of which the thermal state is a special case. One could further generalize the initial state to include 
initial squeezing, etc. For simplicity, we shall not consider here the more general cases. 

As in the previous part we can compute the phase space densities fah for the kink wall, but now for general initial 
densities. Also here we can study the asymptotic late time behavior analytically. Again we are interested in the axial 
vector phase space density fsh, which becomes in the asymptotic late time limit with general initial densities. 



sin(2w+t) 



2w+yP 



(71) 



Thus, compared to Eq. (j60p . there is only an extra factor nh+ — rih- in fsh, which appears similarly in fih and /2/1. 
We can now compute the CP-odd and CP-even phase space densities 



h=± 



^a±h\\fi*±h\'^i 



sin(2a;+t - x (n/j+ - Uh-) 



h=± 



2k 

UJ+ 



-(1 - 2|/3±|-) + ^ — cos(2w+t - A(p) 



X {nh+ - rih-) 



(72) 



where Aip is given in Eq. (j62p . For comparison, in appendix [Cl we compute the fermionic phase space densities fah 
with a general initial state for a thin wall. 
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V. A COMPARISON WITH GRADIENT APPROXIMATION 



In Refs. [4^|4^ (see also Refs. [3a,|33) it was shown that the gradient approximation, when apphed to the evolution 
equations for the Wightman functions, yields a semiclassical force which affects motion of particles in a plasma as a 
(planar) bubble wall of a first order electroweak transition sweeps through the electroweak plasma. The force is of 
the order h, it is proportional to the spin orthogonal to the planar wall, and it has opposite sign for particles and 
antiparticles. Sinc e up to now no analogous analysis has been performed for the time dependent case studied here 
(see however Refs. |59l - l6l| ). we shall present such an analysis here. The Dirac equation for a Wigner transformed ()50|) 
Wightman function =piS'^^(u; v) = {'tp{u)il!{v)) is of the form, 

{7"fco-7-fc + ^7°5t- [77l«(^)+^m,(^)7']exp( - ^UU-o)}*^+~(fc;^) =0, (73) 

where x = (u + v)/2 and A:^ is the Wigner momentum (the conjugate to u — v). One can show that the operator in 

Eq. (j73p commutes with the helicity operator H = k ■ 7*^77^, implying that the Wightman function can be written 
as a sum of helicity diagonal 2 x 2-blocks (|?T|) . With the Ansatz (|?T|) one can construct non-local partial differential 
equations for the densities gah- The real and imaginary parts of these equations yield the constraint equations (CEs) 
and kinetic equations (KEs), respectively. The CEs and KEs can be subsequently solved in a gradient approximation. 
The technical details of these steps are in appendix |F] Here we only present the main results. To first order in 
gradients, the CEs for g^h and g^h do not contain derivatives (here h ~ \): 

(fc2_|m|2-fc2^goh = (74) 

,2 




kt - |m|^ - fc^ - ^^ e\rn\^ )a^k. = 0, (75) 



where the mass has been written as mn + mii = |m|e*^, and k — \\k\\. General solutions to Eqs. ([7Wf5|) are of the 
form, 

ffo/i 5o/t27r(5(fco - |mp - fc^) , ffsfe = g3/t27r(5(^fco ~ ~ ~ ^ ^"^fc"^*^ ) ' ^^^^ 

Equations ([7il - [75|) reveal that at first order in derivatives (a) the vector density g^h does not feel any effect of a 
changing background, while (b) the axial vector density g^h lives on a shifted energy shell given by: 

2; 



\m\'dt 



c^3h = + /i^rf^ ; ^o(t) = v/^-2 + |m(i)P. (77) 

In analogy to the case of a planar wah, in which case the energy shift is, 5ijJ±s = =Fs|™P((^z^)/[2wo-\/^l + \m\'^] , see 
Eq. ([1]) and Ref. [1^, where k± is the momentum perpendicular to the wall and s = ±1 is the corresponding spin 
eigenvalue, Eq. ([77| shows that the axial density g^h lives on a shifted energy shell produced by the zero component 
of a fictitious 'axial-vector field' h{\m\'^dtO)/ {2ku!o)- In this case however, the energy shift is proportional to helicity 
instead, it has identical sign for positive and negative frequency states and, as expected, it is proportional to a time 
derivative of the rate at which the mass argument 6 = Arg[m] changes, which is a (good) measure of CP violation. 
Thus, just like in the case of a planar wall, the time dependent effect is a CP-violating shift, and thus can represent 
a source for baryogenesis. 

In order to determine how this energy shift affects particle densities, we need to consider kinetic equations (to 
second order in gradients). These are derived in Eqs. (jF22HF24| of appendix [F] For goft and g^h the KEs are 

a , dt\m\'^ gah „ ,„„s 
dtgoh + — ^ — ^"~k^ " 

Otgsh + — dk^g^h + h — dk^gsh = . 79 

2ko zkko 

Equation ([75)) teaches us that, as expected, the vector density goh does not feel any force at second order in gradients. 
The only effect that goh feels is a classical 'force', which is of first order in time derivative, and takes account of 
the energy non-conservation in a time dependent background uJo{t). On the other hand, we see from ([7^ that the 
time dependent energy shift effect ([77]) produces as expected a second order semiclassical 'force' term in the kinetic 
equation for g^h- These results are in accordance with what one would expect based on Eqs. (|74H75I) . Just like in 



15 



the planar wall case, there is no second order 5^.2 term; only a term containing single fco derivative occurs in (|79p . 
justifying the name 'force'. In fact, there is no force in ([79| . A better analogy is the Lorentz 4-force F^, where the 
3-Lorcntz force F = e{E + vxB) constitutes the spatial part of F^^, while the 0th component F^ ~ ev ■ E yields the 
rate of energy loss in an electromagnetic field (which of course does not depend on the magnetic field B). Similarly, 
in the above equations we can identify the rate of energy loss as the zcroth component of a 4-force, 

|H!,,Mg!M.ft.,,,,,, ,80, 



where we projected fco — > ±^3/1 on-shell in Eq. (|79p . 

Now, the quantities fah considered in the rest of this paper are simply related to gah via the integral (|52p . In the 
light of Eq. ([76]) we see that the integral ([52]) just projects gah on the positive and negative frequency shells. Unless 
given differently by initial conditions, the positive and negative frequency projections are the same, and this fact does 
not change with time, because the semiclassical force is the same for both frequency shells. This is to be contrasted 
with the planar wall case, in which case the energy shift at the order h has an opposite sign, see Eqs. (HHH). Of course, 
this simple picture is true only in the absence of interactions. When interactions are included, one expects off-shell 
effects in gah, and by performing the integral ([5^ one in general loses information. 

Let us now integrate Eqs. (|78p and (j79p over fcg- Integrating the first equation is easy, and yields a conservation of 
vector phase space density, 

dtkh^^, (81) 

which is consistent with the more general result, which states that /o/i = const in a free theory. Integrating Eq. ([75]) ) 
is more delicate, and yields 

dthh + 1 + h —2 ) hh = . (82) 



Since the expression in the parentheses is dt ln(a;3/i), this equation can be simplified to, 

dtHhh) = -dtHu:^h{t)) , (83) 

and its solution is simply 



/3.(i) = -^/3-., (84) 



where lo- is given in (j3ip and /g^ = fshit — ^ —00) (j70p . This means that, if one starts with f^h = —{kh/uoJ) x {nh+ — 
fih-) (see Eq. ([TO))), the gradient approximation yields, 

kh 

f3h{t) = j-{nh+ - fih-) . (85) 

W3hilj 

This result shows that gradient approximation captures the change in the frequency felt by particles, but does not see 
any quantum effects such as squeezing. Having a cursory look at figures [3HS] shows a striking feature: the large and 
oscillating contribution in the axial vector density is completely missed in gradient approximation. In afterthought, 
this should not come as surprise, since the oscillatory contributions to the densities come from state squeezing, which 
is a genuinely non-adiabatic quantum effect, and thus cannot be captured in a gradient (adiabatic) expansion. The 
question is whether this failure of gradient approximation means that an important effect is missed this way in regards 
to baryogenesis sources. The answer is not so simple as the plots in figures [3HS1 suggest. Note that, when averaged 
over time, the oscillatory contributions disappear, and one is left with a mean effect. This mean effect is captured 
(to a certain extent) in gradient approximation (|85p. Indeed, Eq. (|85p contains a CP-violating contribution, which is 
present during the time transient, and can be extracted from the CP-odd part of Eq. 



E/3'^W = ^#^(l-2nth), (86) 

h=± 

where for simplicity we took an initial thermal state, — hji^ —5- 1 — 2nth, "-th ~ l/[cxp(/3a;_) + 1]. 

In order to compare the CP-odd axial density in the gradient approximation (j86p to the exact results (see figures 
I3HS]), we integrate the phase space densities over the momenta and sum over the helicities, which gives the CP 
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Figure 6: Difference between the finite temperature and 
zero temperature CP-odd axial vector current for the ex- 
act solution (blue, solid) and gradient approximation (red, 
dashed) for a thick wall with 7 = Cl/?""*^. The mass pa- 
rameters are: mi — 0.1/3^^, miR = 0.1/3^^, m2R = ■ 



Figure 7: Difference between the finite temperature and 
zero temperature CP odd axial vector current for the ex- 
act solution (blue, solid) and gradient approximation (red, 
dashed) for a wall with 7 = ■ The mass parameters 
are: mi = 0.1/3~^, m\R = 0.1/3~^, m2R = j3~^ ■ 




Figure 8: Difference between the finite temperature and 
zero temperature CP odd axial vector current for the ex- 
act solution (blue, solid) and gradient approximation (red, 
dashed) for a wall with 7 = 3/3^^. The mass parameters 
are: mi = 0.1/3""', miR = 0.1I3~^, m2R = /3~^. 



odd current ^fi^± jsh- The zero temperature part of the current however diverges as fc — >■ 00. Therefore we only 
compare the finite temperature parts of the integrated fshS, that is, only the part that is Boltzmann suppressed. 
Technically we compute J2h=±[j3h{P) — jzhiP — > 00)], the difference between the CP-odd axial vector current at finite 
temperature and the current at zero temperature. 

In figures [51 [7] and [5] we show the finite part of the CP-odd current for the gradient expansion and for the exact 
solution, for a thick wall with 7 = 0.1/3~^, a wall of thickness 7 = /3~^, and a thin wall with 7 = 3/3~^ respectively. ^ 
The mass parameters are chosen such that there is a large state squeezing. The gradient expansion captures quite 
well the main trend, but misses the oscillations at later times. It is intriguing that already for the wall with 7 = 
the exact solution for the current starts to look quite different from the current in the gradient approximation. For 
the thin wall in Fig. |8]the difference is even more significant. It would be interesting to explore in more detail the 



^ Note that, based on Eqs. (jsjl and Q in the introduction, modes of both the thin and thick wall regime contribute to the integrated 
density, i.e. to the current. The question is therefore, whether the currents in Figs. |6H8]are dominated by modes that satisfy the thick 
wall condition, or by thin wall modes. Roughly speaking, for f}-) <g 1 the current is largely dominated by thick wall modes, whereas for 
/37 3> 1 modes that satisfy the thin wall condition contribute mostly to the current. 
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CP-violating current in the thinner wah regime. Needless to say, in order to make a definite statement about the 
validity of gradient approximation, one needs to perform a more detailed analysis which includes scatterings coming 
from quantum loop effects. 



VI. CONCLUSIONS AND DISCUSSION 



In this work we derive an exact solution of the Dirac equation for fermions with a time dependent mass, generated 
by a scalar Higgs condensate. We assume that the mass has a tanh time dependence, which represents a quite realistic 
model for phase interface (bubble wall) at a first order phase transition in the early Universe setting. Moreover, the 
mass is complex with a phase changing in time, which can be a source for CP violation. We have studied this CP 
violation by looking at the CP-odd part of the axial vector current, since that current biases sphaleron transitions. 

As already emphasized in the introduction the division between "thin wall" and "thick wall" cases depends on the 
relevant momenta, and both cases are present in a typical (say thermal) distribution. For large state squeezing, i.e. 
a large late time Bogoliubov particle production, the axial vector phase space density for a thick wall experiences a 
large oscillatory enhancement and phase shift with respect to the thin wall case. This non-adiabatic behavior cannot 
be captured by the gradient approximation. However, the mean effect for the axial vector current, which is obtained 
from the phase space density after summing the momenta, is described reasonably well in gradient approximation. 
Still, the exact solution for the axial vector current shows that the difference for thinner walls can be quite significant. 
This invites for a more detailed quantitative study of baryogenesis sources in the thin wall regime. 

Extensions to our work can be foreseen. First of all, instead of a time dependent mass, our analysis can be 
generalized to a planar wall case, in which the mass is dependent on one spatial coordinate. This extension is in 
principle straightforward, and has been considered for the CP-conserving case in Ref. [sgI [STj. except that we intend 
to work now with full inclusion of a CP-violating mass. Moreover, whereas we consider one non-interacting fermionic 
species, one could also study multiple flavors. In that case, there is an additional CP-violating source in the flavor 
mixing mass matrix. It would be interesting to see what is the dominant source of CP violation depending on wall 
thickness. Finally, we did not include other effects such as plasma scatterings, and it would be worthwhile to explore 
how these affect the analytical results in this work. 
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Appendix A: Mode functions for a constant mass 

In this appendix we construct the fundamental solutions of Eqs. (|17 p - (|18p for a constant mass. In appendix [B] we 
treat the thin wall case. 

When the mass is constant, the general solution of Eqs. ([T71 - [TS)) (see also Eqs. (fTl fTC)) ) can be written as plane 
waves, 

Lh = A^e^*"* -I- B^e*"* ; Rh = Che-'^' + Dhe'^' 

Lh = Ahc'^' + Bue-'^' ■ Rh = Che'"^' + DhC-''^\ (Al) 

where Ah, Bh, Ch, Dh, Ah, Bh, Ch and Dh are constants. The first order equations ([TTHTB]) tell as that Ch and Dh 
can be expressed in terms of Ah and Bh'. 



„ m* , uj + hk , m* uj + hk , „ m* „ uj — hk „ m* uj — hk „ 

- ^k^'^ = ^^'^ = H V ^^'^ ' = = -^^'^ - - H V ^'"'^ ■ 

(A2) 

Notice that these equations also imply an on-shell condition, ijp = fe^ -I- |mp. Analogous relations hold for the barred 
constants, 



^ m* ^ uj + hk ^ m* uj + hk ^ ^ m* ^ uj ~ hk ^ m* u) — hk 

uj — hk ^ m ^ \m\y (jj — hk ' uj + hk ^ m ^ |m|V^^ + ^fc 



(A3) 
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These constants can be further constrained by imposing the vector current conservation law, which in the absence of 
flavour mixing becomes particularly simple, dti^{x) = 0, or cquivalently j'^{x) = {^'^ {x)il'{x)) = const. In order to fix 
the constant, notice that 

{^Pl{x,t)^Pa{x,t)) = l^[^t(^^i)^^^(^'^i)]^ ^ ^S^p6^ix-x')~F^p{x,t;x',t), (A4) 

where Fai3{x,t;x' ,t) the Hadamard (statistical) Green function, which in the free space vanishes ^. When written in 
momentum space 

^e^'=-^(fc,t) (A5) 

Eq. (jA4p becomes, 

where the expectation value was taken with respect to the vacuum state |f2), in which case F ^ 0. Similarly we have, 

{^l>a{k,t)ii)l{k,t)) = ^Sa^f! + Fc,fs{k,t;t) ^Sc^fs. (A7) 
Taking a trace of (|A6HA7|) one finds, 

Y,{\L,f + \R,f) = iTr[<5„^] = 2 , + \R,\') = l^S^p] = 2 , 

h h 

which implies: 

\U\^ + \Ru\^ = l = \Ln\^ + \Ru\\ (AS) 

where we assumed that the vacuum state normalization is independent of helicity. Moreover, since • £,h' = 5hh', 
Eq. (|A8p agrees with Together with (|A2HA3p the condition (|A8p allows one to completely specify the vacuum 

fermionic mode functions for constant mass (up to an overall phase) , 



L,ik,t) . \/^e— , R,ik,t) . V^M^ 



Notice that these solutions satisfy the mode normalization conditions (|A8[) . Moreover, (when summed over h) they also 
satisfy the (stronger) consistency condition ()11[) . The only remaining conditions to check are the mode orthogonality 
conditions (fT2|) . They imply. 



L*iRh + R*hLh — 

For the solutions above, this implies, 

1 r . 



2uj 



0, 



which can be only satisfied if S[m] = 0. There is no problem when m is time independent. In this case one can 
perform a global rotations on spinors that removes the imaginary part of the mass, 

V' ^ e-^'^'^V = [cos(6') - I sin(6l)7^]V' , V' ^ V^e"'^^' = V;[cos(6l) - i sin(6')75] , 



* This definition of the Hadamard function differs by a constant from the definition used in e.g. Ref. [63 |. which is due to the normal 
ordering of the creation and annihilation operators, assumed in the construction of the density operator in Ref. [63l . 
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where 

tan(20) = — , cos(20) = , sin(20) = ^ , 
TOfl \m\ \m\ 

with which 

— m^p^/jji — m*^l)'ijj]^ = —mii'iijtp — inij^j^^p — — Iml^il! . 

Of course, this rotation is global, and works only if m is time independent. Here we are interested in a time dependent 
problem, and hence the mode functions are not orthogonal as in ^2]). Is this a problem? Not necessarily. What 
is important is that the mode functions span the whole Hilbert space. The most important condition that must be 
satisfied is in the end Eq. ([T3|) . Note also that 

Xhik, t) ■ xh{k, t) = — = -Dh{k, t) ■ Vh{k, t) . (AlO) 

Appendix B: Mode functions for a step fermion mass 

We shall now assume that the mass takes a simple form, such that it exhibits a sudden jump at t = 0: 

m{t) = 7n_e(-t) + m+e(t) (Bl) 

where m_ and m+ are (in general complex) constant masses at negative and positive times, respectively. This is 
what we refer to as the thin wall mass profile. One can easily convince oneself that a constant U{\) rotation of the 
left- and right-handed spinors in ([S]) can remove a jump in the e.g. imaginary part of the mass. After performing 
such a rotation the mass can be written as m± = \Tn± \ exp[i(0-|- -f where x is the relative phase between left- and 
right-handed spinors. One should then solve |m+| sin((/)_|_ -fx) — sin((/)„ -f- x) for x- ^ Therefore, without any 
loss of generality, one can assume that 3[m(i)] = const., or equivalently, 5[m+] = 3[m_]. This will be important for 
decoupling of the equations for u±h defined below (sec also Eq. (US])). 

For the problem at hand with ni{t) given in (jBl[) . wc can make the following Ansdtze for the mode functions (c/. 
Eqs. (|A9|), 

Ln{k,t) = e{~t)L- +9{t)L+ 
Rh{k,t) = e{-t)R^ +e{t)R+ 

Lhik,t) = e{-t)L^ + e{t)L+ 

Rh{k,t) = 9i-t)R- + 9it)R+ , (B2) 
where the solutions for t < are vacuum mode functions derived in (|A9|) . 



uj^hk _t uj^+hk ml - \ uj_ - hk ^ - uj^+hk m*_ ^ 

(B3) 

The solutions for t > can now be written as a linear combination of the normalized positive and negative frequency 



The exact solution gives 



X = arctan 
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solutions 



i , ■. J- LiJ+ — hk „, , ,1 UJ+ + hk , f 



"y 2w+ |m+| 2a;+ |m+| 



where 



2 



and the solutions multiplying fi^ and are the normalized negative frequency solutions of Eqs. ([T7|) - ([T5)) . The four 
Bogoliubov coefficients in (jB2p are determined by the matching conditions. Eqs. (| 171418]) together with the structure 
of the mass term (|B1|) tell us that the mode functions must be continuous at i = 0, which implies the following four 
matching conditions 

Lh{t^Q-)^Lh{t^O+); Rh{t = 0-) = Rh{t ^ 0+) ; = 0-) = = 0+) ; Rh{t = 0-) = Rhit ^ 0+) . 

(B5) 

These conditions imply for the Bogoliubov coefficients in the solutions for Lh and Rh in Eq. (jB2[) : 





- hk 


UJ+ 


+ hk 


UJ^ 


- hk 


UJ+ 


- hk 



UJ+ 


+ hk 


CJ_ 


- hk 


UJ+ 


- hk 




- hk 



By observing that + ~ 1 can easily check that these conditions satisfy the correct normalization 

condition 

To summarize, our simple calculation shows that, as a consequence of a sudden mass change (at t = 0), the number 
of fermions (of each helicity h) produced is 



nh = \P^\ = -. = ^ (B7) 



which is helicity independent and, in general, docs not vanish. One can show that Uh in (jB7P satisfies, < nh {k) < 1 
for arbitrary k, m+ and to_, which is a nontrivial check of the correctness of ()B7p . Note finally that, when fc = 0, 



jm+l |to_ I 



(B8) 



This last result nicely illustrates the dependence of particle production on the complex phase of the mass. Let us 
denote m± = \m±\e^^^. Then 

nh = ^ [1 - cos(0+ - 0- )] = sin^ (^^^^ , (B9) 

which shows that n/j varies between and 1, as it should. There is no particle production when 6^ — 9- = 27rn 
(n € Z), in which case there is no CP violation. On the other hand, when fc = particle production maximizes for 
0+ — 0_ = 7r(27i +1) (n G Z), in which case CP violation is maximal. This is in accordance with the results of finite 
wall thickness (j45l - l47| , shown in figures [TH2] 
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Appendix C: Two point functions for a step fermion mass 



111 this section we consider the evolution of the two point functions with a Heaviside step function mass (jBl[) . Here 
we heavily use Ref. [5^. By making use of (|54p we can find the phase space densities for the thin wall mass profile. 
For t < 0, the phase space densities are given in Eq. ([55]) . Together with Eqs. (jB4[) . Eqs. (jB2[) can be used to compute 
fah for t > 0. 

Instead of calculating the fah by using the explicit expressions for the mode functions, we can also derive the phase 
space densities from the kinetic equations (|53|) . We consider a generalized initial state, for which the fah for t < are 
given by Eq. ([7CT)) . Next we solve the fah for t > from the kinetic equations ([5^ . In order to do so, note first that 
Eqs. (|B4[) and ((54)) imply that the general form of the solutions for i > is. 



fah = "a'l cos{2uj+t) + Pah sin(2w+i) + "fah (a = 0, 1, 2, 3) . 



(CI) 



where Uah, f^ah and ^ah are constants that can be determined from the matching conditions at < = 0. Although 
m experiences a finite jump at t = 0, the structure of the equations of motion (|53p implies that all of fah must be 
continuous at i = 0: 



Iahik,t = Q)^f+^{k,t^Q). 



(C2) 



The first equation in (|53p tells us that the vector particle density /oh cannot depend on time, i.e. that a^h = Pah = 0, 
and we have: 



J Oh 



nh+ + Uh- 



(C3) 



The other three equations in (|53p give nine conditions among the parameters {a^/i, /3,7i, 7^^} {i = 1,2,3). These 
conditions represent a highly degenerate system, such that the following independent conditions remain: 



\aih + mfa2h + hkaah = ; m^/3ih + mjP2h + hkPzh = ; 72/1 



In addition, Pih are related to a.ih as follows: 



hk mt 

Plh = "2/1 H "3/1 



Furthermore, the matching conditions (jC2p result in: 



mi 



:lih ■ 



73/1 



hk 

—lih . 



TO 



R 



hh — "l/i —Oi^h ; P^h 



-otih 



-"2/1 



(C4) 



(C5) 



m 



—2^ {nh+ - rt/i_ ) = "i/i + 71/, 

bJ- 



[nh+ - rih-) = "2/1 + 72/1 

OJ- 



- — {nh+~nh-)^a3h+73h, (C6) 



which, together with Eqs. (jC4HC5p completely specify fih- When the solutions for the aih, Pih and 7^/, are inserted 
into the general form (|Cip . we find 



flh 



fih 



J 3/1 



5R[to_|_to!l] 





r 




mj 


mt fc^ 
f — — 


+ 5R[to+to1 








/ifc 


hk fc2 4 
1 


- 5R[TO+ml] 









cos(2a;-|_<) + 
cos(2aj-|_t) — 



hk 
UI-UI+ 

hk 



(to J — TOj!") 



(to„ — TO^ 



sin(2a;_|-i) 



sin(2a;+i)- 



TO^ fc^ + 5R[to+to'*l] 

x(fi,,+ - Hh-) 
fc^ + 3?[to_|_to' 

x(fi,,+ - fih^) 



9[to+to1] ft,fc fc^+ 5R[to+to.1 
cos(2w+t) H sin(2a;+t) 



(C7) 



(C8) 

■x{nh+-~nh-), (C9) 



It can be checked that the same phase space densities are obtained by inserting the mode functions for t > Q Eqs. (|B4 
in the definitions for fah Eqs. (|B2p . 

The final produced particle number (j56p is 
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For the free vacuum, where nh+ = 1 and fih- = 0, this indeed reduces to the result computed before, Eq. (jB7[) . The 
thermal limit is obtained when n^_|_ — 7%^^ — > 1 — 2nthj 'i-th = 1/ (e^'^+ + 1). 

To get the CP- violating density, which is of relevance for baryogenesis, we finally need to sum the axial density f^h 
over /i = ±, 



h=± 



3h 



2k fc^ + '^[m+m*_] 



2\m. 



|7Ti_ I sinff? 4 



2k 



2k + 3?[TO+m* 



C0s(2w+t) X (n;i+ - Uh-) 



-0-) 



UJ+UJ- 



sin(2a;+i) x {nh+ - Uh-) 



(Cll) 



where in the last equality we used, m± = |m±| exp(i0±). The first current in (ICllI) is CP even, while the latter is 
CP odd. and can be used to source baryogenesis. This latter term is there only when the source of CP violation, 
— 6^ — 6^ does not vanish. The first term in the first equation is the (adiabatic) vacuum contribution (i.e. the 
leading classical term that would survive in the very thick wall limit). 



Appendix D: Mode function normalization 

In this appendix we shall show that the properly normalized mode functions that solve Eqs. ((32|) . and whose indices 
are (|33p . are given by Eq. p4p. Since a±, 6±, c in Eq. (|33p are non-integer, the two independent solutions for x±h are 
the usual ones, and they are of the form 

■"S = ufl^^z^il- zf X2Fi{a±,b±;c;z) 

u'-il = u%z"+'-%l - zf X 2Fi(a± + 1 - c, &± + 1 - c; 2 - c; z) 

= u'il,z''+'-%l-zf+---^-'± X2Fi{l~a±,l-b±;2-c;z), (Dl) 

(1 2) 

where u^^g are the normalisation z-indcpendent constants (which wc shall determine below) and 2^1 (a, b; c; z) denotes 
the Gauss' hypergeometric function, whose series around z = reads, 

, , ab a(a + 1)6(5+ 1) z2 Via + n)T{b + n)T{c) z"" 

2F^{a,b-c;z) = l + —z+ '\ -:- + ■■■+ \, ' y ^ + -- - ■ 

c c(c + 1) 2! I [a)l {b)l [c + n) n\ 

Notice that we have picked the sign of a and /3 in (|30p such that the first (second) fundamental solution in Eqs. (jDl[) 

corresponds to the positive (negative) frequency wave at early times. Notice further that the solutions for v_^'^ ~ ""±1/1 

are obtained simply by flipping the helicity h in u^^if' . The latter form for of the two solutions (|D1[) is useful in that 
the prefactor is in the form a + \ — c = —a = a* , 13 + c — a± — b± = — /? = (3* . It then follows that, in the vicinity of 

t — )- —00 {z e^'''*), u^^l and reduce to the positive and negative frequency solutions, respectively, 

u^^l « 41e^-* = u%e--' , u^ll « 41e-^-* = ^^e^-* {t ^ z ^ e^-*) . (D2) 

In analogy to what we did in Eq. (jA9[) . here we shall take the positive frequency solution at t -> —00, i.e. u±h ~ u'l^l- 
Due to the fact that the vacua at t — >■ —00 and t — -l-oo are not the same (they are related by a Bogoliubov 
transformation), the positive frequency solution becomes a mixture of positive and negative frequency solutions close 
to t ^ +00, as can be implied from the following identity (see e.g. Eq. (9.131.1-2) of Rcf. (63l|). 

2Fi{a±,b±;c;z) = y ^'^ x 2Fj_{a±,b±; a± + b± + 1 - c;! - z) 

r(c- a±)r(c- b±) 

+ ^^^^tM^7I^^(1-^)'""*"'* ^2Fi(c-a±,c-6±;c+l~a±-6±;l-z). (D3) 

r(a±)r(6±) 

Indeed, we have, 

u±h = u'-ll = uihoz^il- z)^ X — vFT A: X 2Fi(a±,6±;a± + 6± + 1 - c; 1 - z) 

r(c- a±)r(c- 6±) 

+ ^±/»o ^^^^T^/°^l',^r ""^ ^"(1 - z)^+--°±-''± X2^^i (c - a±,c- 6±; c+ 1 - a± - 6±; 1 - z), (D4) 



23 



from which wc infer that in the vicinity oft— ^oo(z— >1 



-27^ 



u±h~u±hOTr/ TFT J—^^ +u±ho — — . e + 



r(c-a±)r(c-5±) 



r(a±)r(6±) 



(t -)■ oo, 1 



-27* 



(D5) 



where we made use ofc = a±-|-&±-|-l — c=l-|- 2(3, (3 + c — a± — b±=P + l — c = —(3. 

One can also construct late time positive and negative frequency solutions that solve the differential equation p2p . 
Equation (jD4[) implies that, if (|Dip are solutions, so must be both parts of Eq. (jD4p . such that the two linearly 
independent late time solutions are. 



~(2) 



-,(2) 



±ho^"il ~ z)P+'-^ X 2Fi(a± + 1 - £, 6± + 1 - 5; 2 - 5; 1 - z) 

±(1 - zf+'-' X 2Fi(l - a±, 1 - fe±; 2 - 5; 1 - z) 

13 



_ ,-,(2) a+c-a±-b 



X2^^i(a±,6±;£;l-z) (5=1 + 2/3) 
-a = a* and /3 - 



~ f 1 2) 

where Mj_;jq arc normalisation constants. Now, because a + 1 — c 



asymptotic forms for the mode functions are. 



,-,(1) ^ ,-,(1) 



''±h 



,•-,(.2) 



(t ^ oo, 1 - z e-27t) 



(D6) 

'13 ^ P*, the 

(D7) 



One can check that Eqs. (jD6[) indeed solve Eq. ([5^. so they constitute legitimate linearly independent solutions for 
the mode functions. And moreover, each of the solutions (|D6p can be written as a linear combination of the early 
time solutions (|Dip . as they should. 



Next, we need to properly normalize our mode functions (|D1P and (|D6p . Rather then performing a quantum 
mechanical normalization such as was used in (56| . we shall use the field theoretic normalization (jlllllSp . since it is 
more suitable for baryogenesis applications we have in mind. Since u+h and u-h are related by a first order differential 
equations (|22p . their normalisation constants are not independent. Let us begin by rewriting (j22p as 



\ d mm + TO2i?(l - 2z) 

z(l — z)— ± I 

dz 



27 



u±h 



hk ± imj 
27 ^ 



(D8) 



By making use of the identities, 

c-1 a±+b±-c 

a = : p = 

2 ' 2 

we can recast Eq. (|D8p as 



■ d b±{a± - c) ■ 

^(1 - z)— + z 7-^ - b±z 

dz a± — 1 — o±) 



± 



27 



1 + b± - a± 



± 



irriiB. _ (a± + b± — l)(a± + 6± + 1 - 2c) 
27 ~ 4(a± - 1 - 6±) 



u±ho X 2-Fi(a±,&±;c;z) 



I MzpM X 2i^i(a± - 1,0± + 1; c;z) , (D9) 
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where we used aip = 6± + 1 and bz^ — a± ~ 1. In order to transform the parameters of the hypergeometric function 
on the left hand side such to correspond to those on the right hand side, let us first make use of the following two 
identities (see Eq. (9.137.6) and (9.137.17) in 0), 



dz 



(2^1(0,6; c;z)) 



ab 



X 2^1 (a + 1,6 + l;c+ l;z) 



abz 



X 2-Fi(a + 1,6 + l;c+ l;z) = (c-1) 2-Fi(a, 6; c - 1; z) - 2-Fi(a, 6; c; z) 



upon which Eq. (jD9P reduces to 

(1 - z)(c- 1) X 2Fi(a±,6±;c- l;z) + 



/ b±{a± - c) 
Va± - 1 - 6± 



c+l + (c— 1 — 6±)z X 2Fi{a±,b±;c;z) 



(DIO) 



u±ho (Dll) 



hk ± imj 

I Uzp/iO X 2Fi{a± - 1,6± + l;c;z) 

27 



This can be further transformed by making use of (9.137.17) in |63| . 

(c - 1) X 2i^i(a, 6; c - 1; z) = 6 X 2-Fi(a, 6 + 1; c; z) + (c - 1 - 6) x 2-Fi(a, 6; c; z) 



(D12) 
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into 



(1 - z)b± X 2Fi{a±,b± + l;c;z) + ^^^^^^^ ^ ^ , x 2Fi{a±,b±;c; z) 



a± — 1 — 6-1 



u±ho 

hk ± imi 
27~ 



(D13) 



We need one more transformation |7l| . 

(1 - z){a-b - 1) X 2Fi(a,fo+ l;c;z) (a - c) x 2^1(0- 1, 6 + 1; c; z) + (c - 1 - 6) x 2Fi{a,b; c; z) , (D14) 



with which one gets on both sides of equation (|D13|) a function with identical parameters, implying the following 
relation between the normalisation constants 



u±ho hk ± imj a± — b± — 1 

^ ^ 

u^ho 27 b±{a± - c) 



(D15) 



Several comments are now in order. This expression shows that CP violation is in the relative phase between u+^i 
and u^fi solutions, 



^^g^^^ _ hk ± imi 



(D16) 



which was to be expected, meaning that there is no trace of CP violation in the parameters a±, b± or c of the 
hypergeometric functions. Moreover, 

k"^ + 7nj _ b±{a± - l)(a± - c)(6± + 1 - c) 



472 (a± -b±- 1)2 

such that the ratio (jPlSp can be expressed in terms of the phase dcp and the parameters a± , b± and c. 



u±ho hk±imj \/—b±{a± — c)(a± — 1)(6± + 1 — c) hk±imi uj- ± {miR + ■m2R) 

" , X " ; ^ , X 



b±{a± - c) 



2 4_ \/ cj_ =F {miR + m2B,) ' 



(D17) 



where we made use of {a± — b± — 1)2 = — [±i(a± — 6± — 1)2] = — 4m2^/72 and in the last step. 



b±{a± - c) =F 



[uj_ =F {miR + ■m2R)]m2R 



(a±-l)(5± + l-c) = ± 



[w_ ± {miR + m2R.)]m2R 



Eq. (|D17[) nicely separates the relative CP- violating phase between positive and negative frequency modes and their 
amplitude ratio which is not CP- violating. From Eq. (|D16p it also follows that 9cp = ^'cp+ = ^^'cp-, i-e. that 

giScp+giScp- — 2 

What remains to be done is to perform normalisation of the mode functions. Eqs. (|A8[) and (|2ip imply for the early 
time mode functions (jPip . 

\u+^f + \u^i,\^ = 1 = \v+h\' + \v-h\^ , (D18) 
and analogously for the late time mode functions (|D6p . 

\u+h\^ + \u-h\^ = 1 = \v+h\^ + \i-h\^ ■ (D19) 
Making use of (jPlTp . the first condition in (|D18[) can be written as. 



\u+ho\ 



2Fi(a+, 6+; c; z) x 2^^1(2 - a+, -6+; 2 - c; z) 
b+{a+ - c) 



{a+ - l){b+ + 1 - c) 



X 2-Fi(a_, c; z) x 2i^i(2 - a_, 2 - c; z) 



(D20) 
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where we made use of the fact that a and /3 are purely imaginary, and of aj_ = 2 — a±,bj_ = —b±, c* = 2 — c and 
z* = z is real (the proper sign in front of the second term is a minus). Next, it is convenient to replace a_ and by 



a+ and 



a+ — 1, a_ = + 1) in the second term. The above analysis (see Eqs. (jD9flD15"|) ) implies 

a-b-1 



b{a — c) dz 



and also {a ^ 1 — b, b ^ 1 — a, c ^ 2 — c) 



a-b-l 



(l-a)(c-5 - 1) 



z(l-z)-+(l 



a-b-l 
c-b-1' 



2-Fi(a, b; c; z) — 2-^1(0 — 1, & + 1; c; z) 



iFiil -b,l- a;2- c;z) = 2^^1(2 - a, -fe; 2 - c; z) 



The latter relation can be used to replace the second hypergeomctric function on the first line in (jD20|) . while the 
former to replace the first hypergeomctric function 2^i(a_, 6_; c; z) = 2-^1(1+ ^ 1,^+ + l;c;z) in the second line 
of ((D20)) to obtain. 



1 = \u+ho\^{2Fi{a+,b+;c;z) 



b+{a+-c) 



a+ — b+ — 1 . d , 

[i l-a^)ic-b^-lf ^^-'^Tz+^' 



a_|- — 5+ — 1 



2i^i(l-&+,l-a+;2-c;z) 



a^ — b^ — 1 d . 

(a+-l)(6+ + l-c) [M^^T^"" rf^^ ' ^ 



a+ — b+ — l 



c-b+-l 

2Fi(a+,b+;c-z) x 2Fi(l-a+, l-6+;2-c;z) 



Next we can make use of the Wronskian for the hypergeomctric functions, 

W[2Fi (a, 6; c; z), zi-'=(l - z)^-''-''2Fi(l - a, 1 - 5; 2 - c; z)] ^ (1 - c)z-=(l - z) 



': — a — b— 1 



from which it follows that, 
W[2Fi{a, b- c; z),2Fi(l-a, 2-c; z)] +2^1(0, b; c; z) 



1 — c c — b — a 



2Fi(l-a,l-6;2-c;z)] 



z 1 — z 

When this is inserted into ()D2ip many terms cancel and one ends up with, 

2 (1 - a+)(c - 6+ - 1) _ cj_ + (mi_R + m2fl) 
(a+ — o_|- — 1)(1 — c) 2lj^ 

To summarize, we have found that the normalized early time mode functions (jPlI) arc (c/. Eq. (|D17P ). 



1 - c 



(D21) 
(D22) 

(D23) 
(D24) 



(1) h^- + {miR + 7n2R) Q/1 ^/3 cn / , x 
+h = \ ^ (1 ' ^) ^ 2Fiia+,b+;c;z) 



U-h = u 



_ (1) _ hk - imi - {miR + m2R) 



2uj^ 



X z"(l-z)'^ X 2F1 (a_, 6- ;c;z) 



(D25) 



which are also given in Eqs. (j34p . An analogous procedure yields two other normalised mode functions given in 
Eqs. 



Appendix E: Connecting the kink wall to the thin wall case 

Here we shall show that the early and late time solutions for the mode functions for a tanh wall profile are 
equivalent to those for the thin wall case in the limit 7 — ^ 00. Similarly, we shall demonstrate the agreement between 
the Bogoliubov coefficients a and /3 and the corresponding particle number. 

The early time solutions for the tanh wall (p4| reduce in the limit 7 — > 00 to 



u-n ^ u^ll ^ -^te^ X X e— , (El) 
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where to__r = mm + m2R and m_ = ui-n + imu is the mass for t <Q. The same result can be obtained by solving 
Eqs. (j22p for constant mass m/j = tti^j^ and choosing the positive frequency solution at early time, u±h = u±/joe~*'^~*. 
After normalisation according to jw+^p + = 1 the result (jEip is found. If we solve similarly Eqs. (|22p for the 

negative frequency solution at early time, v±h ~ v±ho&^'^~* , we find that 



V+h = u 



+ (-h) 



V-h 



2w_ 
hk + imj 



— X c 

2tj_ 



(E2) 



In order to compare to the thin wall solutions at early time Eqs. (jA9p we should rotate back to the Lh, Rh basis from 
the u±h basis. By making use of ([2T|) and the solutions (jElHE2p we compute 

)^ — hk + TO_ „ , . + hk + m* 



LI = 



■\/4c<j_(a;_ + TO„fl,) 
— hk — 7Ti_ 



■\/4a;_(cj_ + itl^r) 
+ /i/s — ml 



At first sight these solutions do not seem to be consistent with Eqs 

L- = 



5|. However, they can be rewritten as 



(E3) 



- hk 



2uj- 



LI 



j_ — hk 
2tj_ 



2w_ 



2a;_ 



(E4) 



where the (real) phases are given by 
Ol = arctan 

9i = arctan 



m/ 



cj- — hk + m^R 
—mi 
Lu^ — hk — m^R 



Or = arctan 



= arctan 



-mj 



cj- + hk + m^R 
mi 



hk — m^R 



(E5) 



Thus the early time mode functions Lh arc Rh in the thin wall limit (jE3p only differ from those computed directly 
for the thin wall (jB3p by a common phase factor. A further global U{1) rotation of the Lh,Rh spinor Xh by e~'^^, 
and of the Lh, Rh spinor i^h by e*^'^~^^\ reduces the solutions (|E4p to those in Eq. (jA9p . Here we used that 



|m_ I 



(E6) 



which follows from the identity arctan(a;) — arctan(y) — arctan((a; — y)/(l + xy)). 

Next step is to check the late time solutions. Since the general late time solution for the mode functions is a linear 
combination of the fundamental late time solutions uf^^ and Uj^^ given in Eqs. (p6l - [37|) . we should consider both in 
the thin wall limit 7 — > oo. Eq. psp now becomes 



U+h 



u^h 



a+hU+h + P+hu\h 



2W4 



hk — irai I cj+ — ?7i+ 



2tj, 



The thin wall limit for the Bogoliubov coefficients (jll]) is (see also (|44l) ): 



/ a;+ ± (mi_R + m2fl)] + a;+ T 2m2 ^ 
a;_ [uj+ ± {miR - m2R)] 2uj+ 



I UJ+ [uj- ± {miR + TO2J? )] UJ+ - uj- ± 2m2R 
uj- [uj+ =F {miR - TO2fl)] 2w+ 



(E7) 



(E8) 
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Thus the coefficients are real and they can be shown to satisfy (|42|) . The late time solutions for Lh and can be 
derived via Eq. ([?!]) from the solutions (|E7p . 



I LJ^ — hk + m_(_ 

i/j^ = ot+h — =e 

hk + m* 



+ ft/s + ml , Wj- — hk — m*. t 



As for the solutions for t < 0, we also write the late time solutions in a form more similar to (|B4|) . This gives 



(E9) 



\ w+ — hk 
Q!+/iW — e i e 



2a;4 



V 2w+ y 2w+ 



The corresponding phase factors are 



g(2) 



arctan 



arctan 



TO/ 



CJ+ — hk + m+B. 
—mi 



3(1) 



arctan 



-TO/ 



/l/s — TO+/J 



fl^j^' = arctan 



CJ-)- + /ifc + TO+/J 
TO/ 



— /l/s — TO+/J 



(ElO) 



(Ell) 



We have seen that the early time solutions for Lh,Rh in the thin wall limit, Eqs. (|E4[) . differ from those computed 
directly for the thin wall (|A9|) by a global phase factor. This factor could be removed by rotating the particle spinor 
by a factor of e~*^^. Since the general late time solution matches the early time one, this means that also the late 
time solution should be rotated by the same factor. The resulting phase factor for the a+h and /3+/i solutions should 
match the phase factor present in the solutions (|B4| . Indeed we can show that 



5(1) 

R 

9(2) 
L 

5(2) 
R 



= Arg[a+] 



7L = 



7L = 



7L = 



Arg[a+] 
Arg[/3+] 
Arg[/?+] 



arctan ( ] 



arctan 



—mi 
m+R 



(E12) 



where the and on the righthand side are the ones in Eq. (|B6p . The conclusion we therefore make is the 
following: although the Bogoliubov coefficients computed by taking the thin wall limit of the kink wall solutions (|E8[) 
appear different from those directly computed for the thin wall (|B6p . they in fact only differ by a global phase factor, 
which does not affect the particle number. In the thin wall limit computations, the Bogoliubov coefficients (jESP are 
real, but the early and late time mode functions carry a global phase factor, see Eqs. (jE5|) and (|E11[) . In the direct 
thin wall computations, the (coefficients of) the mode functions do not carry the global phase factor, sec (|B3|) and 
()B4|) . but the global phase is present in the Bogoliubov coefficients ()B6|) . which arc hence complex. In any case, global 
rotations of the (anti)particlc spinors are always allowed, and such rotations simply move the phase factor back and 
forth between Bogoliubov coefficients and mode functions, leading to physically equivalent solutions. 



Appendix F: Deriving the kinetic and constraint equations in gradient approximation 

In this section we derive the kinetic and constraint equations for the densities gah- Our starting point is Eq. ([73 
with the Ansatz ([51]) . After inserting 7^7" = 1 in front of iS^ (fc; x) in Eq. (f73|) . we get, 

{(/ ® /)fc° -p^®{k-cT) + ^iI® I)dt - mR{t){p^ ® I) exp ( - - mi{t){p^ ® /) cxp ( - ^%%o) } 



x{p''gah{k;x))®-{l + hk-a)^0, (Fl) 
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where we made use of the Bloch representation of the Clifford algebra, in which 

7° ^ / ; f ^ ® zcr* ; 7^ -p^ ® I ; H = k- 7^77^ k ■ I ® a . (F2) 
Now, upon multiplying from the left by 

{/, hYj'',~ihf, -7'^} -> {/ (g) /, (K) ha\ p^ ® ha\ p^ (g) /} , 
and taking a trace, we get the following four equations. 



goh - mii[ 



kog2h - ihkgih + ^dtg2h + im_R(i)exp(^ - ^%~^ko'^g3h - mi{t)cxp(^ - 
kogah - hkgoh + T^dtgsh - imii{t)cxp( - ^%l^ko)92h + zm/(i)cxp( - ^%~^ka)9ih 





)92h 


= 


(F3) 




)93h 


= 


(F4) 




)9oh 


= 


(F5) 




]9ih 


= 0. 


(F6) 



Now, hermiticity of ij^S~^ implies reality of the component functions gah, such that the real and imaginary parts of 
Eqs. ()F3HF6| must be separately satisfied. The real parts yield constraint equations (CEs), 



kogoh~hkg3h-mii{t)cos(^-%l^ko^9ih-mi{t)cos(^-%'dko)92h = (F7) 

kogih - mR{t)cos(^^%~^ko^9oh - mj{t)sm(^^%~dko)93h = (F8) 

^052/1 + TOfl;(i)sin(^i^4'^feo^53,i - m/(t)cos(^iVt'^fco^.goh = (F9) 

fco53/i - /ifcffoh - '7iij(t)sin(^iVt'^fco^g2/i + n^/(i)sin(^iVf'^fco^.gi;, = 0, (FIO) 



while the imaginary parts yield the following kinetic equations (KEs 

■l^^ \ /Itr^ 



dtgoh + 2mR{t)sm(^-%l^ko^9ih + 2mi{t)sm(^-%l^ko^g2h = (Fll) 

c^tgi/i + 2/ifcg2h + 2m/?(t)sin(^iVtltfc„^.goh - 2m/(t)cos(^iVtl|:(,^.g3h = (F12) 

9tg2/i - 2/ifcgi;, + 2m/?(i)cos(^iVtl|;„^.g3h + 2m/(i)sin(^iVt'^fc(,^.goh = (F13) 

dtg3h-2mR{t)cos(^^%~dko'^92h + 2rni{t)cos(^^%~dko^9ih = 0. (F14) 



Due to the nonlocal nature of the derivative operators, these equations are hard to solve, and hence not very useful. 
However, when truncated at a finite number of derivatives dt, they reduce to a set of relatively simple equations. This 
derivative truncation, which is a generalization of the quantum mechanical WKB expansion, holds when formally 

h\\dt\\ <t: \ko\ , (F15) 

where we reinserted h to make it explicit that this derivative expansion is in fact an expansion in powers of h. Note that 
no matter how large the norm ||9t||, there always will be modes that satisfy the criterion (jFlSp . Conversely, no matter 
how slow the changes in time are, there always will be modes that break (jF15[) . In some sense, the criterion (|F15[) 
divides a theory into two parts: the semiclassical part where (|F15[) holds and the quantum part, where (|F15[) is broken. 
Of course, a full quantum mechanical kink wall treatment is required for those modes that break condition (jFlSp . 
while a semiclassical treatment should sufhce when (|F15|) is satisfied. When modes are massive on both sides of the 
wall, then the theory has a gap of 2min[|m(t)|], and - at least on-shell - |fco| > min[|m(t)|]. ^ 



When interactions (loops) are included, due to quantum effects the mass gap can decrease, or even completely disappear, so one should 
be careful when making statements concerning applicability of the gradient approximation, even in the case when the tree level mass 
is present on both sides of the 'wall'. For example, in the case of the electroweak phase transition, it is typically the case that the tree 
level mass, and hence the gap, is zero on the symmetric side of the bubble wall. 



29 



While the classical kinetic theory is obtained by keeping the CEs up to 0th order in derivatives and kinetic equations 
to first order in derivatives, in order to get semiclassical equations which contain information on CP violation require 
to maintain first order derivatives in the CEs and second order derivatives in the KEs. Let us now consider the 
constraint equations (|F7HF10[) . We have 



gih 

92h 



mR(t) 



cos(^- dt dkojgoh + sm(^- dt Okojgsh 
I = ; sm[-dtdk„)g3h -\ ; cos - 9t 9fe„ 



(F16) 
(F17) 



Upon inserting these into Eqs. (|F7|) and (jFlOp . and truncating to first order in derivatives, we get 



ko 



-goh 



kogsh 



mudtmi - midtmR 

hk + — dk„ 

2fcn 



iriRdtmi - midtiriR 1 
2 



gsh 
goh 



(F18) 
(F19) 



These two equations can be easily decoupled, such that (again to first order in gradients) we have 



Imp -k^ -hk 



ml 



gsh 




0, 



where we used a shorthand notation. 



ran = \m\ cos{0) ; 



mi = TO 



sin(6') 



(F20) 
(F21) 



mRdtmi - midtmR . 



Eqs. (jF20p and (jF21[) are presented in the main text, see ([7^ and (|75p . Next we consider the kinetic equations to 
second order in gradients. First we treat the kinetic equations for go/i a-nd g^h, hi order to describe CP violation in 
the axial vector current. So, we begin by inserting Eqs. (jF16IIF17|) into Eqs. (|Flip and (jF14p . and we get 



dtgoh + — :^9ko— = 
2 fco 



dtgsh 



dt\m\ 
2/co 



■dkogah 



goh 

ko 



0. 



(F22) 
(F23) 



where we made use of 



mRdtruR + midtmi = -9t|™p 



mRdfmi — midtmR = dt{\m'\^dtO) 



(F24) 



Now, upon making use of 50/1 = {hko/k)gzh plus higher orders (c/. Eq. (jFlSP ) and pulling fco to the left of the 
derivatives, Eq. ()F23P simplifies to. 



otg?,h + -rrr^ — Ofeoffsh + h — — dkogvi = . 



2fco 



2fcfco 



(F25) 



These are presented in the main text in ([78]) and (|79)) . 



Additionally, we can solve the constraint and kinetic equations for gih and g2h- In a similar procedure as 
before, we first take the CEs (|F7p and (jFlOp . and combine them to find 



= (fco - k'^)goh - komRCOs(l-%~^k„) gih + hkmisin(^%~^k„)gih - komicos(l-%l^k„) g2h - hkmRfi\n(\\~^kQ 



g2h 



{kl - k'^)g3h + komism(^-%l^ko^gih - hkmRCOs(^-%~iko'^gih - komRsm(^-%l^ko'^ g2h - hkmicos(^-%l^ko'^ g2h 

(F26) 
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Now, by making use of these equations we can eliminate goh and g^h from the remaining constraint equations (|F8[) 
and (jF9p . After an expansion up to second order in derivatives dko we find, 



kg — k'^ — m|j 



kogih + hk 



2 kl-k^""'' 2 



t.2 iL2 



rriRmR kp mRiTiR fcp , 1^2 
:^ Ci n T";^ TTT H :^ t";^ rrr Ci 



fco 



miTTifj 

f^2 _ f^2 



'^0 y!j,2 _ 1^2 



kog2h + hk 



1 „ mirfii 1 
-f^fcn — Ofc, 



9 h.2 _ u2 '"'■^0 



" kp ~ k'^ 



9ih 

9lh 
92h 



i.2 i.2 



;^2 _ ^2 



^2 

—kog2h + hk 



8 fc-.^-pSfco-i-Z-H^^o^T^p-'^-o 



,92/1 = 



(F27) 



j9k„ + — ^ — dk, 



1 



9 1,2 _ 1,2 ' 



fc2 - /fc2 



.92/1 



-'fco ^,2 _ y^2 



mimfi 

_ f,2 



kogih + hk 



mirhi 1 rnRiTiR 

dk„ + — ^ — dk 



2 fc2 _ ^.2 ' 



" k^ — k^ 



mimR 2 



fco 



mijm/ fco ^2 1^ • 



ko 1.2 jL2 



-mirriRdko 



32/1 



ffi/i = 0. 



(F28) 



As is mentioned above it is only necessary to solve the constraint equations to first order in gradients in order to 
describe CP violation. The reason for expanding to second order will become clear once we discuss the kinetic 
equations. For now we take the constraint equations (jF27|) - (jF28p only to first order in gradients and proceed with a 
description of the decoupling procedure. 

After multiplying Eq. (jF27p with fcp — fc^ — and adding mimu times (|F28p . the zeroth order contribution in g2h 
drops out. Thus, this gives the zeroth order shell for gih, which follows from 



fco 



1,2 jL2 



[(fco - fc^ - m|f)(fco - fc^ - mj) - mjmj^] gih = fco(fco - fc^ - \m\'^)gih . 



Similarly, we can find the zeroth order shell for 32/1 • Although the constraint equations are now decoupled at zeroth 
order, both gui and g2h still contribute at first order. However, at this point we may use the zeroth order relation 
between 52/1 and gn^. From Eqs. (|F27p and (jF28p it can be seen that 52/1 = [n^ii^iR/ {k^ — fc^ — m'f)\gih and 
gih = [mimii/(kQ — fc^ — m|j)]g2/i, respectively Inserting the first in the constraint equation for gih-, and the latter 
in the constraint equation for g2h (both decoupled at zeroth order), one obtains 



kri k 



k^-k'- \m\ 



hk 



rhj\ 



gih = 



hk g2h = . 



gih = gih'^TrS ykl - fc^ - |mp + hk 
52/1 = hh^nS fco - fc^ - |m|^ - hk 



Note that the first order derivatives dkg have been canceled. The solutions for gih and g2h are 

TO/ 

mR 
rriR 
mi 

Thus, like the axial density g^h, also gih and g2h live on shifted energy shells 



(F29) 
(F30) 

(F31) 
(F32) 



uJih = Wo - hk- 



TO/ 



W2/1 = oja + hk- 



ruR 



where ujq is presented in ([77]). Note that the CP-odd part can be shown more explicitly by writing 



(F33) 



LOih = Wo - hk- hktanO- — j — r 

Zluq zwo to 



uJ2h = Wo + hk- hk cot 6- — j — r 

Zluq zwo to 
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CP violation is present due to the changing phase of the mass. 

We continue with the kinetic equations (jF12p and (jF13p . Also here we eliminate gg/i ^nd 93;; in favor of gih and 52/1 
and we expand to second order in gradients. The result is 



dtgih - 
+ 2hk 



miTTiR 
miTTiR 2 1 



^0 a , -a ^0 



"ka 1,2 



^0 

rfiRmi 



fc2 _ 



rriRmi 



i-2 _ 1.2 



2hk 



m/ruf ^2 



-52/j 

1 



-miniR 



mjmi 



Ofeo + mimRpko 



fc2 



^0 



8 fc2-fc2 fco 



-rhRrfiRdko- 



fco 



.91/1 

52?!. 



0*32/1 + 2/lfc--2 r^g2/i + 



' k^ 



niRmR 



2hk 



ruRmi 2 
— ; — a, 



rhimR 



fc2 
1 



fc2 



miTUR 



ko l2 
^0 



fc2 



■ ^2 _ y^2 "fco 



fco .2 



gih 



g2h = 



92h 



fr2 



fcg — fc2 



2/lfc 



mBjnR 2 _ 
8 fc2 



1 



-mRTni-^ -jdfe,, + ■mRmidkoT2 — 

f\jQ K nJQ 







fc2 



+ 



mRTTiR 



1 



A-2 



^fco + 7"^/™/9fco- 



A:2 



fc2-fco ■ 4""""""''°fc2_fc2 



gih 



gih = . 



(F34) 



(F35) 



At this point we can use the constraint equations (jF27[) - (jF28P to replace some of the zeroth order terms in Eqs. 
(|F34[) ~ (|F35[) by terms of higher order in derivatives. The remaining terms in the equations above are the zeroth order 
terms dtgih and dtg2h plus a mix of first and second order terms in gih and g2h- In the same fashion as was done for 
the constraint equations, we can eliminate the first order terms for 52/1 from the kinetic equation for gih by inserting 
the first order solutions for 92/1 from Eq. (|F28| (vice versa for the kinetic equation for g2h)- Then we eliminate the 
remaining second order terms for 172/1 by using the zeroth order terms for g2h from Eq. (jF28P (similarly for the second 
kinetic equation). In the resulting kinetic equations for gih and g2h all double derivatives 32^ have dropped out. The 
result is 



Otgih + — Okogih 



dtg2h 



2ko 
dt\m\'^ 
2ko 



rriRmR 



hkrriRrhi 



dkog2h 



^2 _ ^2 _ 
mjrni 



■2\2 



hkm\dt{mi/mR) 

gih — /,2 — 73 — r3T^fco.9i/i — L) 



hkrajiJiR 



{kl-B-mlf 



g2h 



2/fco(fco - A:2 - mf) 
hkm^dtirriR/mj) 



2 7.2 _2 \^fc-o52/i - . 



2fco(fc^ - k^ -ml) 



(F36) 
(F37) 



The kinetic equations for gih and g2h can be integrated over fco to obtain kinetic equations for the phase space densities 
fih and /2ft,. By making use of the solutions (|F32p one obtains 



dtHfih)^dt\n['^ 

dt H f2h) = 9* In ( ^ 



(F38) 
(F39) 



These equations are easily solved with initial conditions (|67)) . and the solutions are 



r rnR(t) 
Jih{tj 7— (1 - 2nth) 



f2hit) 



wih(t) 

UJ2hit) 



(1 - 2nth) . 



(F40) 
(F41) 



Thus, we have found the phase space densities fih and f2h in the gradient approximation, including CP violating 
effects. It can be shown that fih and /2/1 from Eqs. (|F40p - (|F4ip . together with the solution for f^h (|55|) and foh = 1, 
satisfy the kinetic equations (|53p to the given order in gradient approximation. 



As a final note we mention that the authors of Refs. [6J-|70j have developed a formalism that takes account 
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of the existence of an additional shell at ko = 0, which is permitted by the constraint equations. For this shell the 
constraint equations can be solved up to zeroth order in the gradient expansion, and its solutions carry information 
on quantum coherence. Note that the /cq = shell seems to be outside the validity of the gradient approximation 
(|F15[) . However, a more general condition for the gradient expansion is ^ 1, for which the fco = shell can 

be incorporated. It is interesting that collective phenomena in the plasma can generate a feature that resembles the 
fco = shell (see figures 13 and 14 in Ref. [l^])- However there are also differences; in the weakly coupled regime one 
sees a double peak structure centered around fco = 0, which as coupling increases merge into one broad shell centered 
around fco = 0. This feature is probably not related to the quantum coherent shell at fco = since it is not generated 
by quantum coherence, but instead by collective plasma phenomena. 
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